Inelastic carrier lifetime in graphene 
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We consider hot carrier inelastic scattering due to electron-electron interactions in graphene, as 
functions of carrier energy and density. We calculate the imaginary part of the zero-temperature 
quasiparticle self-energy for doped graphene, utlizing the GqW and random phases approximations. 
Using the full dynamically screened Coulomb interaction, we obtain the inelastic quasiparticle life- 
times and associated mean free paths. The linear dispersion of graphene gives lifetime energy 
dependences that are qualitatively different from those of parabolic-band semiconductors. We also 
get good agreement with data from angle-resolved photoemission spectroscopy experiments. 
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I. INTRODUCTION 

Graphene, a single layer of carbon atoms covalently 
bonded together in a honeycomb structure (as in a mono- 
layer of graphite), was previouly thought to be unstable 
and hence non-existent in a free state. Recently, how- 
ever, Novoselov et al. reported^ that they had succeeded 
in fabricating single graphene sheets. Subsequently, sev- 
eral experimental groups have reported interesting trans- 
port and spectroscopic measurements^^, which has led 
to experimental and theoretical interest in this field that 
is rapidly burgeoning. 

The overlap of the tt z orbitals between neighboring car- 
bon atoms in the graphene plane is accurately described 
by a tight-binding Hamiltonian. Around the K and K' 
points (often called Dirac points) which are at the cor- 
ners of the hexagonal Brillouin zone, the kinetic energy 
term of the Hamiltonian is well-approximated by a two- 
dimensional (2D) Dirac equation for massless particles, 
Ho = —V{){(T x k x + Uyky), where a x and a y are 2x2 Pauli 
spinors and k is the momentum relative to the Dirac 
points (h = 1 throughout this paper). The two compo- 
nents of the spinors correspond to occupancy of the two 
sublattices of the honeycomb structure in a hexagonal 
lattice. This Hq gives a linear energy dispersion relation 
£fc,s = st?o|fc|, where s = +1 (—1) for the conduction (va- 
lence) band. The corresponding density of states (DOS) 
is D(e) = g s g v \e\/(2TTVo), where g s = 2, g v — 2 are the 
spin and valley (i.e., K and K' points) degeneracies, re- 
spectively. Thus, graphene is a gapless semiconductor. 
In intrinsic graphene, the Fermi level lies at the Dirac 
points, but as with other semiconductors it is possible to 
shift the Fermi level either by doping the sample or ap- 
plying an external gate voltage, which introduces 2D free 
carriers (electrons or holes) producing extrinsic graphene 
with gate voltage induced tunable carrier density. The 
Fermi momentum (kp) and the Fermi energy (Ep, rela- 
tive to the Dirac point energy) of graphene are given by 
kp = (Ann/gsgv) 1 / 2 and \Ep\ = v$kp where n is the 2D 
carrier (electron or hole) density. 

Interparticle interactions can significantly affect elec- 



tronic properties, particularly in systems of reduced di- 
mensionality. Moreover, the linear energy dispersion 
around the Dirac points gives condensed matter exper- 
imentalists a unique opportunity to study interaction 
effects on effectively massless particles. In this paper, 
we focus on the effect of electron-electron (e-e) interac- 
tion effects on the imaginary part of quasiparticle self- 
energies, Imp]. From Im[E], we can extract the quasi- 
particle lifetime, which gives information that is relevant 
both to fundamental questions, such as whether or not 
the system is a Fermi liquid, and to possible applications, 
such as the energy dissipation rate of injected carriers in 
a graphene-based device. In particular, an important 
physical quantity of both fundamental and technological 
significance is the hot carrier mean free path, which we 
calculate as a functions of energy, density and in-plane 
dielectric constant. Furthermore, Im[S], being the width 
of the quasiparticle spectral function, is related to mea- 
surements in angle resolved photoemission spectroscopy 
(ARPES). 

The rest of the paper is organized as follows. In sec- 
tion II we discuss the general theory of self-energy of 
graphene. In section III we present our calculated quasi- 
particle damping rate of graphene and compare with the 
damping rate of parabolic 2D systems. In section IV we 
show the detail reuslts of the self-energy, and finally we 
conclude in section V. 



II. THEORY 

We evaluate the self-energy E within the leading-order 
ring-diagram GqW approximation, which is appropriate 
for weak-coupling systems, given by^ 

E s (fe,iw ra ) = -k B T*^2 X! Go, s >(k + q,iu n +iv n ) 

xW{q,iv n )F ss ,{k,k + q). (1) 

Here, T is temperature, s, s' = ±1 are band indices, 
Go is the bare Green's function, u n ,v n are Matsubara 
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fermion and boson frequencies, respectively, W is the 
screened Coulomb interaction, and F ss i{k,k') — i(l + 
ss' cos 9kk>), where Okk' is the angle between fc, fc', arises 
from the overlap of |sfc) and \s'k'). The screened interac- 
tion W(q,iu n ) = V c (q)/e(q,iv n ), where V c (q) = 2ne 2 /Kq 
is the bare Coulomb potential (k =background dielectric 
constant), and e{q,iv n ) is the 2D dynamical dielectric 
function. In the random phase approximation e(q, iv n ) = 
1 — V c (q)Tlo(q, iv n )i where the irreducible polarizability 
IIq is approximated by the bare bubble diagram*^ , which 
gives the familiar Lindhard expression [with a modifica- 
tion to include the form factor F ss > (fc, fc')]. 

The self-energy approximation described by Eq. 
should be an excellent approximation for graphene 
since graphene is inherently a weak-coupling (or "high- 
density" in parabolic-band systems) 2D system^. The 
measure of the ratio of the potential to the kinetic en- 
ergy in graphene is r s — e 2 /(KVo) « 2.2/ k, where k 
is the effective dielectric constant of the graphene layer 
and the media surrounding it. (Note that, unlike regular 
parabolic-band systems, r s in graphene does not depend 
on density.) Typically in graphene, r s < 1, implying that 
the system is in the weak-coupling regime. 

After the standard procedure of analytical continua- 
tion from iu> n to u> + i0 + , the retarded self-energy can 
be separated into the exchange and correlation parts 
£ s et (k,w) = £° x ( k ) + £™ r (k,w)£ The exchange part is 
given by 

Ef(k)=-£] n F (£ k+q y)K(q)i^(k,k + q), (2) 

s'q 

where np is the Fermi function, and Zj-^ s — e^.s — is the 
electron energy relative to the chemical potential /x. This 
term in graphene is discussed in Ref. d. At typical exper- 
imental densities, the Fermi temperature Tp = \Ep\/ks 
in graphene is very high compared to the temperature of 
the sample {e.g., for n = 10 13 cm- 2 , Tp w 3.1 x 10 3 K, 
and Tp cx yn)- Therefore, it is an excellent approxima- 
tion to set the temperature T = 0, which we do in the 
rest of this paper. This implies that np(^ s ) = #(£ks) 
(where 9 is Heaviside unit step function), and pL = Ep. 

The correlation part, £g° r (k,a->), is defined to be the 
part of T,l ct (k,Lu) not included in £f (k). In the G W 
approximation, the T,° s or (k, uj) can be decomposed into 
the line and pole contributions, S cor = S lino + £P° le io 
where 
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s'q 



do/ y e (q)fW(k,k + q) 
2tt £k+q,s' - UJ - ioj' 



e(q,iu)') 



(3) 



s'q 

x K(q)^ss'(k,k + q) 



-^k+q,s')] 



e(<7,Ck+ q ,s' - oj) 



The E ox and E lme are completely real, the latter because 
e(q, iuo) is real. Thus, Im[£P ole (k, uj)\ gives the total con- 
tribution to the imaginary part of the self-energy, i.e., 



x V c {q)lm 



£(<7,6e+<},s' - W) 



F ss ,(k,k + q). 



(5) 



The inverse quasiparticle lifetime (or, equivalently, the 
scattering rate) r s (fe) of state \sk) is obtained by setting 
the frequency in imaginary part of the self-energy to the 
on-shell (bare quasiparticle) energy £ s fc, i.e., 



r s (fc) = 2Im[£f(fc,£ fes )]. 



(6) 



(4) 



(The factor 2 comes from the squaring of the wavefunc- 
tion to obtain the occupation probability.) The GqW 
self-energy approximation used here is equivalent to the 
Born approximation for the scattering rate. Note that 
the integrand of Eq. ([5]) is non-zero only when Im[e] 
or Re[e] = 0. These correspond to scattering off single- 
particle excitations and plasmons, respectively. 



III. QUASI-PARTICLE SCATTERING RATE 

The self-energies and quasiparticle lifetimes of 
graphene and conventional parabolic-band semiconduc- 
tors differ considerably. These differences can be ex- 
plained with the help of Fig. 1, which shows the single- 
particle excitation (SPE) and injected-electron energy 
loss (IEEL) continua and the plasmon dispersion for a 
direct gapless 2D parabolic band semiconductor and for 
graphene. The intersections of the IEEL continua with 
the SPE continua and the plasmon dispersion lines indi- 
cate allowed decay processes via e-e interactions. In both 
doped parabolic-band semiconductors and graphene, an 
injected electron will decay via single-particle intraband 
excitations of electrons in the conduction band. In 
parabolic band semiconductors, an electron injected with 
sufficient kinetic energy can also decay via plasmon emis- 
sions and intcrband SPE (also known as "impact ion- 
ization"). On the other hand, as shown in Fig. 1(b), 
electrons injected into doped graphene cannot decay via 
plasmon emission, and the region in q—u> where decay 
via interband SPE is allowed (along the straight line seg- 
ment between uj = v^kp and vo(k — kp), for k > 2fcp) is 
of measure zero. In fact, within the Born approximation 
the quasiparticle lifetime of graphene due to e-e inter- 
actions at T — is infinite. (Multiparticle excitations, 
which are excluded in the approximations used here, will 
give a quasiparticle a finite lifetime^, but this is a rela- 
tively small effect in graphene.) 

In doped graphene, the only independent parameters 
relevant for Born appxoimation quasiparticle scattering 
rates at T = are the Fermi energy relative to the Dirac 
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FIG. 1: (Color online) The single-particle excitations (in- 
traband and interband) and injected-particle energy-loss 
(hatched region) continua, and plasmon dispersion (thick line) 
for (a) gapless parabolic-band semiconductor with equal hole 
and electron masses and (b) graphene, at T = with Fermi 
energy in the the conduction band. For gapped semiconduc- 
tors, the interband continua are shifted up by the energy gap. 



point Ep = vokp and the dimcnsionlcss coupling con- 
stant r s = e 2 /(Kv ). The scattering rate, which has 
units of energy, must therefore be proportional to Ep, 
and must be a function only of e/Ep = k/kp and r s . 
Fig. 2(a) shows the Born approximation T — quasipar- 
ticle lifetime 1/r = T due to e-e interactions as a function 
of energy £(k) = efc — Ep. Since the speed of the quasi- 
particles close the the Dirac points is approximately a 
constant vq « 10 8 cm/s, the inelastic mean free path I 
is obtained by = I'o'KC)- In Fig. 2(b), we provide 
the corresponding I, which shows that at n = 10 13 cm -2 
a hot electron injected with an energy of 1 eV above Ep 
has an I due to e-e interactions that is on the order of 
20 nm. This will have implications for designing any hot 
electron transistor type graphene devices. In particular, 
because Klein tunneling 1 — in graphene creates problems 
in the standard gate-potential switching method in tran- 
sistors, in its place could be a switch based on modulating 
the electron energy £ a regime where \dtjd£\ is large. 

As with doped parabolic-band 2D 
semiconductor ;) 13 ' 14 ' 15 , in graphene T(fc) = l/r(k) cx 
(k — kp) 2 \ \og(\k — kp\)\ for k w kp due to scatter- 
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FIG. 2: (a) Inelastic quasiparticle lifetime/scattering rate 
(1/r = r) in graphene due to dynamically screened e-e inter- 
actions, as a function of energy at T — for different densities, 
within the Born approximation, (b) The corresponding quasi- 
particle mean free path for n = 10 13 cm -2 (corresponding to 
Ef ~ 0.4 eV). (c) Equivalent scattering rate for a parabolic 
band semiconductor (without interband processes), for com- 
parison. 



ing phase-space restrictions. Further away from kp, 
however, the dependence of T on k in graphene and 
in parabolic-band semiconductors are markedly and 
qualitatively different. To wit, in parabolic band semi- 
conductors plasmon emissio n 14 ' 16 and interband collision 
thresholds 1 ^ cause discontinuities in the T(k), as shown 
in Fig. 2(c), but the graphene T(fc) is a smooth function 
because both plasmon emission and interband processes 
are absent. 



4 



IV. SELF ENERGY 

In order to see the effects of the plasmons and in- 
terband SPE in graphene in Imp rct ], one must look 
oif-shell, i.e., uj ^ £ fe>s . The off-shell Imp ret ] is not 
merely of academic interest; it is needed to interpret 
data from ARPES. The spectra of the ARPES electrons 
ejected from graphene give the electronic spectral func- 
tion, from which one can infer E ret (fc, u>}^. Physically, 
in the GqW approximation, the off-shell 2Imp rct (k, uj)] 
gives the Born approximation decay rate of the quasipar- 
ticle in state k if it had kinetic energy uj + Ef rather than 

In Fig. 3 we show Imp ret (fc, oj)] for k = and k = kp 
as a function of uj. Within the GqW approximation, the 
contributions to the off-shell Imp rct (fc, uj)] can be visu- 
alized as the intersection of the SPE continuum and plas- 
mon line in Fig. [I] with the vertically displaced IEEL 19 . 
At k — there are two contributions to Imp], the intra- 
band and interband SPEs. For low energies (|cj| < E F ) 
only intraband SPE contributes to Imp]. Its contribu- 
tion reaches a maximum around Fermi energy, then de- 
creases gradually with increasing energy, as is the case 
with a parabolic-band semiconductor— where it is the 
only decay channel for the quasiparticle (assuming uj < 
band gap energy). But in graphene there is a new de- 
cay channel of the quasiparticle, the interband SPE. Due 
to the phase space restrictions the interband SPE does 
not contribute to the self energy at low energies, but 
at higher energies (oj > Ep) its contribution increases 
sharply, overwhelming the SPE mtra contribution. The 
gpgmtcr con t r ibution then increases almost linearly with 
oj, with the same slope as for intrinsic graphene^. Plas- 
mons do not contribute to Imp rct (fc = 0, oj)] for uj > 0. 

At k = kp, not only do plasmons contribute to Imp], 
in the low-energy (oj < 2Ep) regime, they actually dom- 
inate over the SPE contributions, as can be seen in 
Fig. 3(b). In contrast, in parabolic-band 2DEGs the plas- 
mon and SPE contributions to Imp rct (fc^, oj)] go as uj 2 
and uj 2 lnw, respectively 1 -, and hence both contributions 
are roughly equal in magnitude. 

Fig. 4 shows the imaginary part of the quasiparticle 
self-energies for the conduction band, Im[E™'], and the 
valence band, Imp™*], of graphene for several different 
wavevectors. For k > kp, the Imp™*] shows a sharp 
peak associated with the plasmon emission threshold. In 
the high energy regime, the dominant contribution to 
Imp] comes from the interband SPE, which gives rise 
to an Imp] is linear in uj for all wavevectors. Note 
that within the GqW approximation, for a given r s in 
graphene, the E, oj and k scale with Ep, Ep and kp, 
respectively; i.e., for fixed r s , the function E(fc, oj) is uni- 
versal, where E = Y*/Ep, k — k/kp and to — uj/ Ep. 

Finally, Fig. 5 shows the calculated Im[E rct ] at fixed 
k = kp as a function of energy uj for various densi- 
ties. No fitting parameters were used in this calcula- 
tion. These results compare favorably with the recent 
data from ARPES experiments by Bostwick, et al£ from 
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FIG. 3: (Color online) The imaginary part the retarded self- 
energy at T = (thick line) for (a) k = and (b) k = k,F as 
a function of energy. The dot-dashed and dashed lines in (a) 
are the SPE intra and SPE intor contributions, respectively. The 
dot-dashed and the dashed lines in (b) are the total SPE (in- 
traband and interband SPE) and the plasmon contributions, 
respectively. 



which they extract Imp ret (fcp, oj)] for different densities. 
(The overall scale is different because we have assumed 
a Si02 substrate used by some other groups, whereas 
Ref. d's samples were on SiC, which has a different n.) 
Bostwick, et al. invoked plasmons, SPE and phonon ef- 
fects in interpreting their data. We find that includ- 
ing just the plasmon and SPE effects, we get reasonable 
agreement with their data, except for features near 0.2 eV 
which probably can be explained by calculations that in- 
clude electron-phonon interactions. 

Before concluding, we discuss some of the approxi- 
mations that lead to Eq. (JXJ) . We ignore on-site (Hub- 
bard) e-e interactions because, at zero magnetic field, 
this interaction is irrelevant in the renormalization group 
sensed. Our calculation does include both intra- and 
inter-band scattering processes, since our expression for 
E s (k, iu> n ), Eq. |T]), contains a pscudospin sum s' = ±1. 
Finally, there is an issue about the coupling of the bands 
due to interactions; i.e., because of non-diagonal terms 
E ss / with s ^ s' , the elementary excitations are a super- 
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FIG. 4: (Color online) The imaginary part of the self en- 
ergy as a function of energy for wavevectors k — 0, 0.5, 1.0, 
1.5/cf for (a) the conduction band and (b) the valence band 
of electron-doped graphene at T — 0. 
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FIG. 5: (Color online) Scattering rate as a function of energy 
for different densities, n — 1.2, 3.0, 5.6xl0 13 cm -2 (for clarity, 
successive lines are shifted upward by 0.05 A -1 ), in the form 
of Fig. 3 of Bostwick et aZ.— . 



position of states from the conduction and the valence 
band. Dyson's equation in this case is G _1 = Gq 1 — S, 
where the quantities are 2x2 matrices, which, when 
written out in full, is 



-E_+(k,w) 



-S+_(k,w) 
w-£_(A)-E__(k,w) 



(7) 



The excitation energies are given by the poles of G, 
or equivalently, the zeroes of the determinant of G . 
Our calculation is equivalent to ignoring the off-diagonal 
terms S |_ in the determinant of G . The off- 
diagonal terms £±=p are each of order r s , and therefore 



the product S-| E 



•"s which can be ignored in the 



weak-coupling regime. We also ignore intervalley scatter- 
ing because the Born approximation Coulomb scattering 
rate in 2D is cx q~ 2 , implying that small-g intravalley 
scattering processes dominate over the large-g interval- 
ley ones. 



CONCLUSION 



To conclude, we have calculated the electron-electron 
interaction induced hot electron inelastic scattering in 
graphene, finding a number of intriguing and significant 
differences with the corresponding 2D parabolic disper- 
sion systems. Our infinite ring-diagram GqW appraco- 
mation should be an excellent quantitative approxima- 
tion for graphene since graphene is a low-r s (i.e., weak- 
coupling) system. Our T = results should remain for 
T <C Tp, which is typically the case in experiments. We 
obtain good agreement with recent ARPES data with- 
out invoking any phonon effects. ARPES data has re- 
cently created a controversy with respect to the role of 
plasmons^. Our detailed calculation using a realistic 
model generally agrees with the interpretation of the data 
in Ref. y. The calculated inelastic scattering length I as 
a function of energy £ is by itself is of interest in the 
context of ballistic hot electron transistor applications 
of graphene, where the performance limitation is inelas- 
tic scattering. Furthermore, the Klein tunneling effect 
in graphene makes it difficult to switch transistors by 
a gate-voltage induced depletion (as in a conventional 
MOSFET). In its place, one could imagine a graphene- 
based transistor in which switching is achieved by mod- 
ulation of the injection carrier energy in the regime of £ 
where \d£(£)/d£\ is large. 

This work is supported by US-ONR. 
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